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As  fuel  cells  become  more  prominent,  new  manufacturing  and  production  methods  are  needed  to  enable 
increased  volumes  with  high  quality.  One  necessary  component  of  this  industrial  growth  will  be  the 
accurate  measurement  of  the  variability  of  a  wide  range  of  material  properties  during  the  manufacturing 
process.  In  this  study,  a  method  to  detect  defects  in  fuel  cell  catalyst  layers  is  investigated  through 
experiment  and  mathematical  simulation.  The  method  uses  infrared  thermography  and  direct-current 
electronic-excitation  methods  to  detect  variations  in  platinum-containing  catalyst-layer  thickness  with 
high  spatial  and  temporal  resolution.  Data  analysis,  operating-condition  impacts,  and  detection  limits  are 
explored,  showing  the  measurement  of  defects  on  the  millimeter  length  scale.  Overall,  the  experimental 
and  modeling  results  demonstrate  great  potential  of  this  technique  as  a  nondestructive  method  to 
measure  defects  that  is  amenable  to  use  on  roll-to-roll  manufacturing  lines. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Materials  and  systems  advances  continue  to  improve  the 
performance  and  durability  of  polymer  electrolyte  membrane  fuel 
cell  (PEMFC)  systems,  but  cost  continues  to  be  a  barrier  to  wider 
adoption  [1],  However,  PEMFC  systems  are  currently  being 
deployed  in  some  early  markets,  and  as  this  commercialization 
expands,  there  must  be  a  transition  from  low-volume,  largely 
manual  manufacturing  methods,  to  high-volume,  continuous,  and 
automated  processes.  In  parallel,  quality-control  techniques  must 
be  developed  that  are  sufficiently  sensitive,  rapid,  noncontact,  and 
nondestructive,  and  that  can  be  implemented  in  continuous  or 
automated  production  lines  to  inspect  the  entire  material.  Scalable 
quality-control  methods  are  especially  needed  for  membrane 
electrode  assembly  (MEA)  components,  including  membranes, 
electrodes,  and  gas  diffusion  media  (GDM),  which  are  often 
inspected  visually  before  use  in  a  stack  [2],  Of  particular  impor¬ 
tance,  given  the  high  cost  of  Pt  and  other  noble  metal  constituents, 
is  the  detection  of  variations  in  catalyst  layers.  Catalyst  layer 
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thickness  variations,  such  as  a  spot  with  excessive  or  limited 
loading,  can  yield  performance  losses  beyond  those  expected  from 
the  Pt-catalyst  variation  alone,  especially  at  the  local  area  near  the 
defect  [3-5],  This  is  most  likely  due  to  uneven  compression  of  the 
GDM  at  the  location  of  the  defect,  which  increases  resistance  and 
leads  to  water  management  and  other  issues  [5], 

Today  many  MEA  components  are  fabricated  using  roll-to-roll 
processes  [2],  wherein  a  sheet  of  material  that  is  continuously 
formed  by  coating,  casting,  or  other  methods  is  conveyed  through 
required  process  steps  by  a  roller  system,  and  wherein  the  final 
product  is  wound  into  a  continuous  roll  (see  [6,7]  for  details  on 
general  roll-to-roll  processing).  A  large  existing  industry  provides 
in-line  inspection  devices  for  typical  roll-to-roll  products  such  as 
adhesive  tape  and  paper  towels.  These  devices  typically  concentrate 
on  optical  inspection  of  visible  defects  or  point  measurements  of 
coating  thickness.  Although  these  measurements  can  be  extremely 
useful,  they  cannot  necessarily  be  applied  to  PEMFC  MEA  compo¬ 
nents,  because  (1 )  not  all  defects  are  visible,  particularly  relative  to 
the  typically  black,  nonreflecting  electrode  coatings  and  GDMs;  and 
(2)  point  measurements  provide  only  statistical  data  -  they  gener¬ 
ally  cannot  locate  discrete  defects,  even  if  scanned  across  the  sheet 
of  material  being  fabricated,  as  is  often  done  in  the  industry.  Given 
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that  discrete  defects  can  act  as  failure-initiation  points  for  MEAs  in 
operation  [8],  an  inspection  technique  must  be  developed  that 
measures  every  location  on  the  sheet  material  such  that  defects  in 
MEA  components  can  be  removed  prior  to  assembly  into  complete 
cells.  This  100%  inspection  requirement  is  necessary  to  avoid 
premature  cell  failures,  which  often  can  cascade  and  cause  complete 
stack  failure  and  necessitate  expensive  stack  replacement. 

For  catalyst  layer  defects,  X-ray  fluorescence  (XRF)  [9,10]  has 
recently  been  applied  to  the  in-line  measurement  of  catalyst 
loading  [11],  Unfortunately,  this  technique  cannot  at  this  time 
provide  100%  inspection  for  catalyst  layer  defects.  Acquisition  times 
provided  by  commercially  available  systems  to  provide  for  quan¬ 
titative  analysis  of  Pt  loading  are  on  the  order  of  10  ms  [12],  At  these 
acquisition  times,  assuming  a  roll  width  of  18  in.  and  a  web  speed  of 
30  ft  per  minute,  data  are  gathered  only  every  1.42  ft  in  the  direc¬ 
tion  of  substrate  motion.  Therefore,  this  technique  can  provide  only 
a  statistical  metric  of  Pt. 

Infrared  (IR)  thermography  may  enable  complete,  rapid  (with 
response  times  on  the  order  of  1  s  for  large  areas),  noncontact,  and 
nondestructive  detection  of  defects.  In  fact,  several  patents 
mention  using  IR  thermography  to  detect  various  PEMFC  parame¬ 
ters,  including  MEA  and  membrane  defects  [13—17],  However,  little 
public  information  is  available  about  the  ability  of  specific  excita¬ 
tion  techniques  to  enable  the  detection  of  specific  (size  and  type) 
defects  in  PEMFC  materials.  Although  IR  thermography  techniques 
may  not  be  element  specific  (as  is  XRF),  this  work  shows  that  IR 
techniques  can  be  used  to  identify  catalyst  layer  thickness  varia¬ 
tions  rapidly,  at  a  relatively  low  cost,  and  with  100%  inspection. 

In  this  paper,  a  DC-based  excitation  technique  is  examined 
experimentally  and  theoretically  for  detecting  variations  in  catalyst 
layer  thickness  via  IR  thermography  (henceforth,  IR/DC).  As  an  initial 
study  in  the  use  of  IR/DC  detection  for  catalyst  layer  defects,  this  work 
focuses  on  detecting  square  spots  of  different  thicknesses  in 
a  surrounding  layer  of  nominal  thickness.  Assuming  that  Pt  is 
uniformly  distributed  on  the  carbon  support  in  the  catalyst  layer,  this 
method  allows  for  the  quantification  of  the  areal  Pt  loading  using  the 
known  Pt/C  ratio.  Section  2  describes  the  experiment  and  theory  for 
the  technique.  Section  3  presents  the  results  of  IR  measurements 
with  defined  defects  that  are  used  for  model  validation.  Section  4 
uses  experiments  and  modeling  to  examine  the  impact  of  different 
variables  of  the  technique,  including  excitation  time  and  voltage,  and 
defect  size  and  aspect  ratio,  as  well  as  fuel  cell  testing  data  to  verify 
that  the  technique  does  not  damage  the  electrode  material. 

2.  Experimental 

Catalyst  layer  samples  were  prepared  from  a  catalyst  ink  con¬ 
taining  Ion  Power  Inc.  DE2020  solution  (20  wt%  1100  EW  Nafion), 
n-propanol,  water,  and  TKI<  Inc.  TEC10E50E  46%  Pt/Ketjen  carbon 
catalyst.  The  ionomer  to  carbon  mass  ratio  was  8:10.  The  total  solids 
content  in  the  ink  was  6.8%.  All  components  were  hand  mixed  with 
a  glass  stir  rod  in  a  glass  vial,  sonicated  for  4  h  (in  an  ice  bath),  hand 
mixed,  transferred  to  a  20-mL  syringe,  and  horn  sonicated  for  30  s. 
The  catalyst  ink  was  then  sprayed  onto  a  150-pm  thick,  5-cm  x  5- 
cm  ethylene  tetrafluoroethylene  (ETFE)  decal  using  a  home-built 
ultrasonic  spray  system  utilizing  Sono-Tek’s  impact  line  of  spray 
heads  [18],  The  spray  system  includes  a  syringe  pump,  a  120-kHz 
micro-bore  spray  head  and  power  system,  and  a  programmable 
x—y  stage.  Details  of  optimal  conditions  based  on  spray  material  can 
be  found  in  [18],  The  flow  rates  of  the  catalyst  ink  during  spraying 
and  the  spray  head  conditions  were  optimized  for  the  system  to 
achieve  a  Pt  loading  of  about  0.04  mg  cm-2  in  a  single  coat.  The 
spray  process  was  performed  on  a  heated  stage  at  80  °C.  No  addi¬ 
tional  time  for  drying  was  needed  between  coats.  Defects  were 
made  by  applying  a  removable  polymeric  mask  over  a  square 


region  (typically  1  cm  x  1  cm)  in  the  center  of  the  25-cm2  layer.  The 
mask  was  removed  during  the  deposition  process  to  achieve  the 
desired  (lower)  defect  area  thickness  (loading).  No  residual  adhe¬ 
sive  remained  after  the  mask  was  removed. 

Fig.  1  shows  schematics  of  the  experimental  setup  and  the 
experiment.  Unless  stated  otherwise,  a  21 -V  DC  bias  was  applied  to 
the  setup  using  a  Keithly  2400  source.  For  alternative  voltages  a  BK 
Precision  XLN  10014  Power  Supply  was  used.  Catalyst  layer  samples 
on  the  ETFE  decals  were  fit  between  two  brass  clamps,  which  were 
clamped  with  two  bolts  torqued  to  10  cm-kgf.  To  detect  the 
temperature  variation,  an  IR  camera  (Jenoptik  Vario-Cam  HiRes) 
was  used.  The  camera  has  a  640  x  480  pixels  detector  with  accu¬ 
racy  under  ±1.5  K.  A  uniform  emissivity  of  0.95  was  used  and  was 
determined  empirically  by  comparing  a  piece  of  3  M  brand  masking 
tape,  which  has  a  known  emissivity  of  0.95  [19],  at  the  same 
ambient  conditions  as  the  catalyst  layer. 

For  cell  testing,  the  defect-containing  decals  were  hot-pressed 
onto  Nafion  212  membranes  at  130  °C  under  120  lbf-cm2actiVe  area 
pressure  for  5  min  and  assembled  into  cells  using  SGL25BC  GDLs. 
Performance  curves  conducted  after  break-in  were  run  with 
hydrogen  and  oxygen  utilizations  of  83%  and  50%  (1.2  and  2  stoich), 
respectively.  Forward  and  then  backward  scans  of  the  current  from 
0  to  996  mA  cm  2  were  executed  in  6  increments  with  15-min 
current  holds  at  each  point  along  the  curve.  The  last  5  min  of 
each  data  point  were  averaged  for  reporting. 

3.  Modeling 

The  development  and  analysis  of  the  diagnostic  technique  was 
aided  by  mathematical  modeling.  For  the  computational  models,  the 
governing  heat  conduction  and  current  conservation  equations  are: 

Picw(g[)  -v-(fciVT)  =  Q.  (1) 

-V-(<J,-V$)  =  0  (2) 

respectively,  where  pt,  cPt,-,  <ju  and  k,  are  the  density,  specific  heat 
capacity,  electronic  conductivity,  and  thermal  conductivity  of 
species  i,  respectively,  T  is  the  absolute  temperature,  t  is  time,  and  <t> 
is  the  electronic  potential.  The  Joule  heating  rate  per  unit  volume  is 
computed  from  the  dot  product  of  the  electric  field  and  current 


density,  i, 

0i  (3) 

where  the  subscripts  i  in  Eqs.  (1) — (3)  can  be  catalyst  layer,  ETFE 
decal,  or  catalyst  layer  defect. 

The  boundary  conditions  for  the  model  are: 

T  =  Tr0om  and  V  =  Vapplied  for  dQe  x  =  0  (4) 

T  =  Troom  and  V  =  0  for  3 Qex  =  w  (5) 

n-(-/qVT)  =  hside(T-Troom)  for  3 Qey  =  0  and  d  (6) 

n-J  =  0  for  3i2ey  =  0  and  d  (7) 

n-(-/<ETFEVT)  =  hbottom(T  -  Troom)  for  dQez  =  0  (8) 

n-(-kCLVT)  =  h[op(T-Troom)  for dQe z  =  h  (9) 


where  w,  d,  h  are  the  width,  depth,  and  height  of  the  computational 
domain,  respectively,  and  h^e,  htop,  hbottom  are  the  convective  heat 
transfer  coefficients  for  the  side,  top,  and  bottom  walls,  respec¬ 
tively.  The  parameter  values  used  in  the  numerical  simulations  are 
listed  in  Table  1,  where  the  convection  heat  transfer  coefficients  are 
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fitting  parameters.  In  general,  the  convection  heat  transfer  coeffi¬ 
cient  for  air  is  5-25  W  m-2  K-1  for  free  convection  and 
10-200  W  m-2  K-1  for  forced  convection.  Therefore,  the  convec¬ 
tion  heat  transfer  coefficients  used  in  this  study  indicate  that  the 
heat  transfer  from  the  catalyst  layer  was  combined  free  and  forced 
convection,  which  is  reasonable  considering  that  no  significant 
measures  were  taken  to  isolate  the  experimental  apparatus  from 
the  environment. 

For  the  modeling  domain,  3D  and  2D  domains  were  used;  3D 
provided  higher  resolution  and  2D  provided  faster  computation. 
Fig.  2  shows  a  schematic  depiction  of  both  computational  domains 
(same  as  Fig.  1  with  the  brass  clamps  removed).  The  thermal  effect  of 
the  brass  clamps  was  accounted  for  using  an  isothermal  boundary 
condition,  and  their  effect  on  the  numerical  temerature  profile  was 
found  to  be  negligible  because  the  defect  is  located  at  the  center  of 
the  sample.  For  the  2D  model,  the  catalyst  layer  and  ETFE  decal  layer 
were  combined  to  form  a  composite  layer  (as  shown  in  the  middle 
part  of  Fig.  2)  and  then  converted  into  a  2D  domain  (as  shown  in  the 
bottom  part  of  Fig.  2).  An  averaging  scheme  based  on  volume 
weighted  properties  was  used  to  ensure  that  the  results  between  the 
two  models  were  consistent,  as  discussed  in  Appendix  A.  Although 
the  2D  model  is  significantly  faster  than  the  3D  model,  it  provides 
the  temperature  distribution  of  the  catalyst  layer  top  surface  only.  If 
the  temperature  variation  along  the  thickness  of  the  catalyst  layer  is 
important,  the  3D  model  is  required. 

The  governing  equations  were  solved  using  COMSOL  Multi¬ 
physics.  The  3D  computational  domain  was  discretized  into 
a  tetrahedral  mesh,  whereas  the  2D  computational  domain  was 
discretized  into  a  triangular  mesh,  and  adaptive  mesh  refinement 
was  used  to  obtain  better  convergence.  For  both  cases,  second- 
order  Lagrangian  elements  were  used  with  the  built-in  direct 
solver  MUMPS.  The  solutions  were  considered  a  converged  solution 
when  the  preset  tolerance  value  was  below  10  fi  and  were  mesh 
independent. 

Table  1 

Parameters  used  in  the  numerical  simulations  for  the  base  case. 

Parameter  Value  Unit  Source 


[26] 

[27] 

[27] 

[28,29] 

[28,29] 

[28,29] 

Fitting  parameter 
Fitting  parameter 
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Fig.  3.  False  color  thermal  response  of  a  pristine  sample  (100%  loading),  100%  reduction  (bare  spot),  75%  reduction,  and  50%  reduction  in  catalyst  layer  thickness  in  the  defect  area. 

4.  Results  and  discussion 

4.1.  DC  excitation  and  model  validation 

Fig.  3  shows  1R  thermographs  of  samples  with  1-cm2  square 
defects  of  different  thicknesses  (loading),  excited  under  potentio- 
static  control  at  21  V  for  5  s.  The  current  flows  from  left  to  right  in 
each  image.  Thinner  areas  had  lower  overall  temperatures  than  the 
thicker  (standard)  areas.  The  catalyst  layer  can  be  idealized  as 
a  mesh  of  resistors  in  which  the  current  flow  is  governed  by  Ohm’s 
law.  In  this  scheme,  the  thinner  areas  lead  to  higher  local  in-plane 
resistances,  resulting  in  lower  current  flow  and  temperature 
change.  Because  of  the  higher  local  resistances  in  the  defect  area, 
current  preferentially  flows  around  it,  leading  to  a  local  tempera¬ 
ture  maximum  just  outside  -  in  this  case  above  and  below  -  and 
a  temperature  drop  inside  the  defect  area.  The  temperature  varia¬ 
tions  of  the  samples  are  statistically  significant,  i.e.,  outside  the  90% 
confidence  values  of  pristine  samples. 

Fig.  4  shows  temperature  line  scans  corresponding  to  the 
vertical  lines  in  Fig.  3.  These  scans  highlight  quantitative  differences 
that  can  be  exploited  to  employ  the  IR  thermography  detection 
technique.  The  data  allow  three  observations: 


Under  potentiostatic  control  the  overall  temperature  of  any 
defect-containing  sample  is  less  than  that  for  a  pristine  sample. 
The  temperature  differential  near  the  edge  of  the  defect 
increases  with  the  severity  of  the  defect  because  more  current 
flows  around  rather  than  through  the  more  electrically 
resistive  defect  area.  This  can  be  used  to  identify  a  defect,  as  it 
depends  on  many  process  variables,  including  defect  size. 
Various  data  analysis  methods  may  lead  to  the  identification  of 
a  defect.  For  example,  one  can  examine  the  absolute  value  of  the 
temperature  differential,  or  one  can  look  at  standard  deviations 
along  a  temperature  line  scan.  Fig.  5  illustrates  an  example  for 


one  possible  data  analysis  method,  which  assesses  the 
maximum  temperature  differential  over  each  0.2-cm  increment 
along  a  line  scan  (hereafter  referred  to  as  AT/Ax).  Any  value  that 
is  statistically  relevant,  i.e.,  higher  than  the  99%  confidence 
limits  around  the  pristine  data,  detects  the  defect.  The  figure 
demonstrates  that  for  a  5-s,  21 -V  excitation,  the  detection  size 
limits  for  a  10%  and  a  50%  thickness  (loading)  reduction  are 
1  cm2  and  0.0625  cm2,  respectively.  We  have  previously  shown 
that  no  detectable  performance  losses  are  seen  for  1-cm2  zero¬ 
loading  defects  25-cm2  CCM  [  11  ].  Our  data  suggests  that  we  will 
be  able  to  detect  all  catalyst  loading  variation  defects  that  would 
result  in  detectable  performance  losses. 


1.01  0.02 
Distance  aero: 


;  electrode  /  m 


Fig.  4.  Temperature  profiles  across  a  25-cm2  catalyst  layer  decal  after  5  s  of  21  -V  DC 
excitation.  All  samples  contain  a  single  1-cm  x  1-cm  defect  in  the  center  of  the  sample 
with  100%  (bare  spot),  75%,  or  50%  reduction  of  the  total  catalyst  layer  thickness  in  the 
defect  area.  Error  bars  represent  90%  confidence  intervals  for  excitations  on  multiple 
pristine  samples. 
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Fig.  5.  Maximum  temperature  differential  over  a  0.2-cm  spot  size  (with  respect  to 
position)  across  defect-containing  and  pristine  samples  after  5  s  of  21 -V  excitation. 
Error  bars  represent  99%  confidence  values  about  the  pristine  samples.  The  x-axis 
indicates  the  length  of  one  side  of  the  square  defect  and  the  percentage  of  catalyst  in 
the  defect  area  for  each  sample. 

Fig.  6  shows  the  comparison  between  model  and  experiment  in 
terms  of  temperature  profiles  (Fig.  4)  for  two  excitation  times.  The 
error  percentages  between  the  simulation  and  data  are  within  a  few 
percent  or  less  near  the  defect  region,  thus  validating  the  model.  The 
numerical  results  do  deviate  from  the  experimental  data  at  the  edge 
of  the  samples,  and  this  discrepancy  can  be  explained  by  the 
assumption  of  the  Pt-loading  profile.  An  XRF  line  scan  for  a  sample 
with  a  25%  thickness  reduction  in  the  defect  area  shows  the  Pt  to  be 
nonuniform,  whereas  the  model  assumes  a  step  change  at  the  defect 
and  slight  linear  decreases  away  from  the  defect.  The  assumption  of 
identical  catalyst  thickness  profiles  (defect  and  nondefect  areas) 
along  the  x-direction  is  not  necessarily  valid  and  may  explain  the 
deviation  between  the  data  and  the  model. 


4.2.  Detection  method  and  defect  and  variable  exploration 

As  noted,  the  response  time  for  implementing  a  diagnostic  tool 
as  a  quality  control  technique  for  fuel  cell-related  roll  goods  is 
critical.  The  response  time  of  the  IR/DC  method  was  explored 
experimentally  and  with  the  computer  simulation.  A  1  °C  temper¬ 
ature  variation  was  chosen  as  the  metric  for  a  successful  defect 
detection  and  serves  as  an  upper  detection  time  bound.  This  metric 
is  within  camera  resolution  and  is  supported  when  comparing  to 
the  AT/Ax  and  a  standard-deviation  criterion,  as  it  is  outside  the 
averages  and  confidence  intervals  for  both.  Fig.  7  shows  experi¬ 
mental  and  simulated  data  for  the  time  required  to  reach  a  1  °C 
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Fig.  6.  Comparison  of  3D  simulation  results  with  experimental  data  of  the  tempera¬ 
ture  profiles  across  the  catalyst  layer  after  5  s  and  10  s  of  21 -V  DC  excitation. 
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Fig.  7.  Time  for  1  °C  temperature  difference  to  be  detected  between  the  nondefect  and 
defect  areas.  The  solid  line  shows  the  prediction  based  on  a  2D  simulation  for  a  1-cm2 
defect  with  21 -V  DC  excitation. 


temperature  difference,  using  the  AT/Ax  data  analysis  method,  to 
be  detected  between  the  pristine  catalyst  layer  and  the  defect  area. 
Excellent  agreement  between  experimental  and  simulation  data  is 
shown.  At  the  given  excitation  conditions,  a  bare  spot  can  be 
detected  in  0.2  s,  whereas  a  20%  reduction  in  catalyst  layer  thick¬ 
ness  in  the  defect  area  requires  approximately  2  s  for  detection. 

The  standard  deviation  of  the  temperature  along  a  line  scan  can 
also  be  used  as  an  analysis  method.  The  experimental  error  for 
pristine  samples  is  shown  as  a  black  line  in  Fig.  8.  Values  above  this 
line  can  be  considered  statistically  significant  for  positive  defect 
detection.  The  model  was  used  to  investigate  this  method  for 
different  defect  sizes  and  thicknesses.  The  results  demonstrate  that 
a  20%  reduction  in  catalyst  layer  thickness  can  be  detected  in  less 
than  5  s  for  a  1-cm2  defect  size.  However,  for  defects  with  smaller 
areas,  only  larger  catalyst  layer  thickness  variations  will  produce 
a  detectable  response.  In  fact,  Fig.  8  shows  that  the  standard 
deviation  grows  in  a  nonlinear  and  almost  exponential  fashion  as 
the  thickness  in  the  defect  area  decreases.  Using  the  AT/Ax  instead 
of  the  standard  deviation  data  analysis  method  allows  for  a  faster 
detection  of  a  smaller  defect  (see  Fig.  7).  The  comparison  between 
Fig.  7  and  Fig.  8  highlights  differences  between  the  AT/Ax  and 
standard  deviation  methods,  and  emphasizes  the  importance  of  the 
method  used  for  analyzing  the  raw  temperature  data.  The  actual 
data  analysis  method  used  on  a  fuel  cell  roll  good  may  further 
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Fig.  8.  Simulated  standard  deviation  in  temperature  profile  with  defect  loading  after 
5  s  of  21-V  DC  excitation  for  different  defect  areas.  The  solid  line  indicates  the 
experimental  standard  deviation  in  studies  with  pristine  samples. 
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Fig.  9.  Measured  detection  time  for  1  °C  temperature  rise  as  a  function  of  applied 
potential  for  a  0.25-cm2  defect  with  50%  thickness  reduction.  Open  symbols  depict 
corresponding  simulation  results. 

depend  on  the  process  control  strategy  and  will  be  refined  with 
respect  to  the  kinds  and  sizes  of  defects  that  are  expected  on  the 
manufacturing  line.  In  actual  operation,  the  time  scale  to  perform 
these  mathematical  operations  is  expected  to  be  such  that  both 
could  be  used. 

Another  critical  process  variable  is  the  DC  excitation  voltage. 
This  was  explored  through  simulation  and  experiment,  and  results 
are  shown  in  Fig.  9,  where  a  catalyst  layer  with  a  0.25-cm2,  50% 
thickness  reduction  was  excited  with  different  voltages.  In  this 
figure,  the  time  for  detection  is  defined  as  the  time  at  which  AT/Ax 
along  a  line  normal  to  both  the  defect  edge  and  current  flow  was 
greater  than  1  °C,  which  was  outside  the  confidence  levels  deter¬ 
mined  for  the  pristine  samples;  this  value  is  averaged  over  time 
(0.18  s)  to  account  for  noise.  The  data  indicate  that  excitation 
voltages  above  20  V  result  in  detection  times  less  than  1  s,  and  that 
above  30  V  there  is  no  additional  gain  in  detection  speed.  The 
results  shown  in  Fig.  9  indicate  that  the  IR/DC  method  presented 
here  is  a  valid  candidate  for  rapid  detection  of  catalyst  layer 
thickness  variations  in  fuel  cell  roll  goods. 

The  validated  model  is  run  to  determine  the  detection  limits  in 
terms  of  defect  size  and  shape.  Fig.  10  shows  the  temperature  rise  at 
the  transition  between  the  nominal  catalyst  layer  and  an  area  of 
50%  thickness  reduction  after  1  s  and  3  s  for  a  21 -V  DC  excitation  as 
a  function  of  defect  size.  The  smallest  defect  size  that  can  be 
detected  within  1  s  using  the  AT/Ax  data  analysis  method  is 
~  0.017  cm2.  However,  by  increasing  the  excitation  time  to  3  s, 
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Fig.  10.  Simulated  temperature  rise  with  time  for  a  21-V  DC  excitation  as  function  of 
50%  thickness  (loading)  square  defect. 


Fig.  11.  (a)  Temperature  profiles  across  the  electrode  and  (b)  temperature  rise  between 
the  defect  and  nondefect  areas  and  standard  deviation  of  the  temperature  profile  as 
a  function  of  aspect  ratio  after  1  s  of  21-V  DC  excitation  for  0.5-cm  wide  defects.  The 
inset  in  (a)  shows  the  definition  of  the  aspect  ratio. 

a  defect  smaller  than  0.01  cm2  can  be  detected  with  the  same 
criterion.  The  results  also  demonstrate  a  nonlinear  dependence  of 
the  temperature  rise  with  defect  size,  and  indicate  that  the  method 
is  applicable  over  a  wide  range  of  defect  sizes. 

The  effect  of  aspect  ratio  (Width/Height  of  defect)  for  a  0.5-cm 
wide  defect  was  also  explored  using  numerical  simulation. 
Fig.  11a  shows  temperature  profiles  across  the  electrode  along  the 
dashed  line,  as  indicated  in  the  inset  for  different  aspect  ratios.  The 
temperature  profiles  change  significantly  with  aspect  ratios,  but 
the  maximum-to-minimum  variation  in  the  temperature  profile 
increases  and  then  decreases  with  increased  aspect  ratio.  This 
effect  can  be  seen  in  Fig.  lib,  where  the  solid  line  shows  the 
maximum  temperature  variation  between  the  defect  and  non¬ 
defect  areas.  Conversely,  the  standard  deviation  of  the  tempera¬ 
ture  profile  monotonically  decreases  with  aspect  ratio,  as 
indicated  by  the  dashed  line.  Thus,  using  both  methods  at  the 
same  time  can  enable  not  only  detection  of  a  defect,  but  also 
(potentially)  its  shape. 

4.3.  Effects  of  DC  excitation  on  fuel  cell  performance 

A  key  attribute  of  any  quality  control  technique  is  that  it  must  be 
nondestructive  to  the  sample.  Thus,  it  is  critical  to  know  whether 
the  high  voltage  applied  during  the  measurement  affects  the 
performance  and  properties  of  the  catalyst  layers.  For  example,  it  is 
known  that  carbon  corrosion  (oxidation  of  carbon  in  the  catalyst 
layer)  occurs  in  an  operating  PEMFC,  which  can  lead  to 
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Fig.  12.  Hydrogen/air  polarization  curves  for  pristine,  10-s,  and  100-h  DC-excited 
MEAs.  Inlet  conditions:  100%/50%  relative  humidity,  80  °C,  1  atm  backpressure,  1.2/ 
2.0  stoich.  Error  bars  represent  90%  confidence  intervals  based  on  repeat,  pristine 
MEAs. 


performance  decrease  due  to  loss  of  electrochemically  active 
surface  area  (ECSA),  when  the  cell  potential  is  held  above  0.8  V  for 
periods  of  time  [20—24],  To  test  the  impact  of  the  applied  21 -V 
floating  potential,  a  series  of  pristine  and  DC-excited  samples 
were  tested  using  polarization  curves  [25].  Fig.  12  shows  hydrogen/ 
air  polarization  curves  for  pristine  and  DC-excited  catalyst-coated 
membranes  (cathode  side).  Before  being  hot  pressed  into  an  MEA, 
the  catalyst  layers  were  subjected  to  21  V  in  air  at  ambient 
temperature  and  humidity  using  excitation  durations  of  either  10  s 
(the  maximum  time  used  in  the  quality  control  technique)  or  100  h 
(to  accelerate  any  effects).  As  illustrated  in  the  figure,  neither  the 
10-s  hold  nor  the  100-h  hold  of  the  DC  bias  appeared  to  affect  the 
MEA  performance  significantly.  Table  2  gives  the  ECSA  of  the  same 
set  of  samples,  and  shows  that  the  ECSA  of  the  voltage-held 
samples  was  not  diminished. 

Both  the  mass  and  thickness  of  the  catalyst-coated  decals 
excited  for  100  h  were  measured  before  and  after  excitation.  No 
changes  were  observed.  Carbon  corrosion  can  often  lead  to  evolu¬ 
tion  of  carbon  oxides  and  can  also  lead  to  formation  of  glassy 
carbon,  which  is  much  denser  than  the  porous  carbon  of  a  pristine 
catalyst  [24],  Additionally,  the  in-plane  resistance  of  the  decals  did 
not  change  over  the  100-h  experiment.  Thus,  no  evidence  of  carbon 
corrosion  or  performance  loss  was  found.  When  applied  within  the 
operating  conditions  reported  here,  the  technique  is  noninterfering 
and  nondestructive. 


Table  2 

ECSA  for  pristine  and  DC-excited  MEAs. 
Pristine  samples 


Average 
90%  confidence 


ECSA  (mV1  Pt) 
63.1 


63.4 

62.1 

63.2 

13 


DC-excited  samples 


ECSA  (mV1  Pt) 
68.8 

69.3 

74.3 


5.  Conclusions 


In  this  paper,  a  nondestructive,  rapid  technique  to  measure 
catalyst  layer  thickness  variations,  which  can  be  correlated  to  Pt- 
loading  variations,  was  introduced.  The  technique  involves 
applying  a  DC  excitation  across  a  catalyst  layer  and  measuring  the 
thermal  response  using  IR  thermography.  Experiments  and 
modeling  demonstrate  that  this  technique  can  monitor  small  vari¬ 
ations  in  thickness  with  temporal  resolutions  around  1  s  or  less.  This 
technique  was  examined  to  understand  the  impact  of  defect  size  and 
shape,  detection  time,  and  applied  voltage.  Results  show  that  higher 
voltages  (<20V)  and  longer  times  (~3  s)  improve  the  defect  size 
detection  resolution  to  0.01 -cm2  areas.  In  addition,  various  possible 
control  and  detection  strategies  were  discussed.  We  conclude  that 
this  technique,  in  combination  with  electrified  rolling  electrodes, 
holds  great  promise  for  use  as  an  on-line  quality  control  method. 

Acknowledgments 

This  work  was  funded  by  the  Assistant  Secretary  for  Energy 
Efficiency  and  Renewable  Energy,  Office  of  Fuel  Cell  Technologies,  of 
the  U.S.  Department  of  Energy  under  contract  numbers  DE-AC36- 
08-G028308  (NREL)  and  DE-AC02-05CH11231  (LBNL).  P.K.D.  also 
gratefully  acknowledges  financial  support  through  a  Natural 
Sciences  and  Engineering  Research  Council  of  Canada  postdoctoral 
fellowship. 

Appendix  A 


To  decrease  the  computational  cost  of  the  simulations,  a  2D 
effective  property  model  was  used  instead  of  the  full  3D  one.  In  this 
2D  model,  the  properties  can  be  defined  as: 

„eff  CL  +  fETFE  O' ETFE  (A1, 

aCCL  =  - Z — T7 -  (A.  1) 

tCL  +  ‘ETFE 


0.eff  =  fdefectOCL  +  ^ETFEOETFE  (A2) 

defeCt  ^defect  +  fETFE 

where  Oq~l  is  the  effective  electrical  conductivity  of  the  composite 
catalyst  layer  (combined  layer  of  ETFE  and  catalyst  layer),  o^efect  *s 
the  effective  electrical  conductivity  of  the  defect  area,  tc l,  tETFE. 
tdefect  are  the  thicknesses  of  the  catalyst  layer,  ETFE  layer,  and 
defect,  respectively,  and  ccl,  oetfe  are  the  electrical  conductivity  of 
the  catalyst  layer  and  ETFE  layer,  respectively.  The  effective  thermal 
conductivities  can  also  be  estimated  using  a  similar  approach, 
whereas  the  effective  specific  heat  capacities  for  the  composite 
layer  and  the  defect  region  are  defined  as: 


CaCCL 


frlLPeCp,CL  +  tETFEPETFEcp,ETFE 
f CL  PE  +  fETFEPETFE 


(A.3) 


_  EdefectPoVCL  +  fETFEPETFE^p,ETFE  (A4) 

p, defect  fdefectPCL  +  fETFEPETFE 

respectively.  The  2D  model  uses  exactly  the  same  set  of  governing 
equations  as  mentioned  above,  where  the  subscripts  i  represents 
the  effective  value  of  either  the  composite  layer  or  the  defect  region 
as  described  in  Eqs.  (1 )  and  (2).  However,  the  source  term  in  Eq.  (3) 
is  modified  to  account  for  the  convective  cooling  from  the  top  and 
bottom  surfaces.  The  source  terms  for  the  2D  model  are  written  as: 

Qccl  =  _i.V0  -  ht°P(T  ~  +  frbottomtr  ~  7-room)  (A_5) 

hIL  +  tETFE 


90%  confidence 


Qdrfert  =  -i-yg  -  htop(T - Troom)  +  /lbotlpm(7' - Troc 

tdefect  +  t£TFE 


